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ABSTRACT 

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 9 1 , we have performed aperture 
synthesis 345 GHz continuum and CO(3-2) observations with the Submillimeter Array resolution) and 

high-resolution imaging of polarized intensity at the K^-band using the HiCIAO instrument on the Subaru Telescope 
(Cy.'25 resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 and 
170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner 
holes. The model fitting analysis of the spectral energy distribution reveals an H 2 mass of 2.4 x 10“^ Mq in the cold 
(T < 30 K) outer part at 65 AU < r < 170 AU by assuming a canonical gas-to-dust mass ratio of 100, although 
a small amount (>3 x 10“^ Mq) of hot (T ~ 180 K) dust possibly remains inside the inner hole of the disk. The 
structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not 
directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3-2) emission with a 
velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating 
gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after 
the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks 
studied thus far. 
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1. INTRODUCTION 


(Evans et al. 2009). Thus, studying the structure and evolution 
of these protoplanetary disks provides information about how 
and when planets form. Initially, the disks are optically thick, 
producing broadband infrared (IR) emission well in excess 


Nearly all newly formed stars are surrounded by disks of 
gas and dust, which provide the building blocks of planets 
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of the stellar photosphere. The star+disk system has a rather 
flat IR spectral energy distribution (SED) slope (e.g., Lada & 
Wilking 1984; Lada 1987). By ~5 Myr, nearly all stars lack 
evidence of warm circumstellar dust and accretion onto the star 
(e.g., Ledele et al. 2010; Currie & Sicilia- Aguilar 2011). This 
implies that most protoplanetary disks around solar/sub solar- 
mass stars have disappeared by this time and that the disk 
material has been accreted/dispersed because of processes such 
as photoevaporation or the disk material has been incorporated 
into planets. 

Transition disks bridge the gap between these endpoints 
of disk evolution because their excess emission at some IR 
wavelengths is intermediate between those of an optically thick 
disk and a star lacking a disk (e.g., Strom et al. 1989). Although 
the term “transition disk” includes a diverse set of morphologies, 
a major subset of these objects include those with near/mid 
IR deficits (relative to an optically thick disk) but optically 
thick emission at other, typically longer wavelengths, which 
is indicative of inner holes or gaps (e.g., Calvet et al. 2002). 
IR SED modeling implies that the hole or gap sizes for these 
objects typically range between "^1 and "^50 AU (e.g., Espaillat 
et al. 2010; Merm et al. 2010) and enclose the planet-forming 
region in our own solar system. SED modeling has identified 
numerous other transition disks in nearby star forming regions 
(e.g., Cieza et al. 2010, 2012b; Merm et al. 2010; Currie & 
Sicilia- Aguilar 2011; Romero et al. 2012). Because formation 
of gas giant planets creates holes or gaps in disks, transition 
disks with SEDs consistent with these features may be excellent 
laboratories for studying planet formation. 

Although SED modeling can provide indirect evidence of 
an inner hole or gap, high-resolution observations with long 
baseline (sub-)millimeter interferometers and large aperture 
IR telescopes provide direct evidence of their existence (e.g.. 
Brown et al. 2009; Hughes et al. 2009; Andrews et al. 2011; 
Isella et al. 2012; Cieza et al. 2012a; Mathews et al. 2012; 
Thalmann et al. 2010; Hashimoto et al. 2012; Mayama et al. 
2012). Since millimeter and micron- sized dust grains contribute 
to the majority of millimeter and NIR emissions, it is important 
to compare the detailed spatial distributions of these tracers for 
understanding the physics of the transition disks. 

Sz 91 is an M0.5 young star surrounded by a transition disk 
located in the Lupus III molecular cloud {d = 200 pc; Comeron 
2008). The stellar position, (o'j 2 ooo, ^J 20 oo), is (16^07"^1P.6, 
— 39°03'47('2). The stellar mass and age have been estimated 
to be 0.49 Mq and 5 Myr, respectively (Hughes et al. 1994). 
The Ha line width at 10% of the peak has been measured to be 
283 km s“\ suggesting accretion with 10“^^ Mq yr~^ (Romero 
et al. 2012). Sz 91 is classified as a wide binary with a separation 
of 9" (Melo 2003), corresponding to 1800 AU. However, such 
a wide binary system is unlikely to disturb the disk evolution 
of a host star (Kraus et al. 2012; Harris et al. 2012). Moreover, 
the difference between their proper motions implies that these 
sources are not co-mo ving (Roeser et al. 2010). We therefore 
treat the Sz 91 disk as a circumstellar disk around a single star 
system throughout this paper. 

The SED of Sz 91 is characteristic of the transition disk. The 
following remarkable features of this source can be identified 
in the SED: no significant IR excess with a spectral index from 
Ks to 24 /xm of —2 (i.e., class III in the IR categorization), 
presence of a large dip of approximately 20 /xm, and very steep 
flux density rising between 24 and 70 /xm (Evans et al. 2009; 
Romero et al. 2012). Romero et al. (2012) categorized Sz 91 
as a giant planet-forming disk on the basis of the following 


features: very steep increase in flux density at 24 /xm, which 
indicates a sharp edge at the inner radius of the disk; a clear sign 
of mass accretion onto the central star; a relatively massive disk 
(-5 X 10-3 Mq). 

Although most transition disks with inner holes have flux 
deficits restricted to 1-10 /xm and optically thick emission at 
longer wavelengths, Sz 91 differs from them because its flux 
deficit extends to significantly longer wavelengths (24 /xm), 
and it has strong far-IR to millimeter emission. The total 
flux density of the 870 /xm (345 GHz) continuum emission 
of Sz 91 has been estimated to be 34.5 zb 2.9 mJy (Romero 
et al. 2012). Although the 1.3 mm continuum emission had 
not been detected previously (Nuernberger et al. 1997), the 
recent wide field imaging survey with the AzTEC receiver on 
the Atacama Submillimeter Telescope Experiment (ASTE) has 
clearly detected a 1.1 mm flux density of 27.2 zb 6.0 mJy 
(R. Kawabe et al. 2014, in preparation). These submillimeter 
flux densities imply that the Sz 91’s disk is substantially massive 
with respect to those of other class III sources in nearby 
star forming regions (Andrews & Williams 2005, 2007); the 
flux densities roughly correspond to 2-5 X 10 3 Mq if the 
canonical gas-to-dust mass ratio of 100 is assumed. Such an 
object having significant submillimeter flux density, but no NIR 
excess appears to be at the transition phase from class II to 
III and is quite rare (~5% Andrews & Williams 2005, 2007). 
Although a recent high-resolution imaging survey for transition 
disk objects (Andrews et al. 2011) has successfully revealed a 
large inner hole in the disk up to a radius of ~70 AU, there 
is selection bias for the most massive transition disks. Because 
the submillimeter flux density of Sz 91 is the lowest among the 
transition disk objects studied thus far, considering the distance 
to the source (Eigure 10 in Andrews et al. 201 1), Sz 91 is a crucial 
target for investigating the disk evolution from the perspective 
of a variety of submillimeter flux densities, i.e., the disk mass. 
However, the disk structure has not been resolved thus far, and 
high-resolution imaging observations are urgently required. 

In this paper, we present the first subarcsecond resolution 
images toward Sz 91 at submillimeter and near-infrared wave- 
lengths. In Section 2, observational parameters and calibration 
for the obtained data are described. In Section 3, high-resolution 
images are shown. In Section 4, the detailed disk structure is dis- 
cussed on the basis of model fitting to the SED and the CO(3-2) 
line profile. The unique nature of Sz 91 is determined from a 
comparison with other transition disks. Einally in Section 5, the 
key results of this study are summarized. 

2. OBSERVATIONS AND DATA REDUCTION 

2.7. High-resolution Imaging of the 345 GHz 
Continuum and CO( 3-2 ) Emission 

Interferometric observations of the 345 GHz continuum and 
CO(3-2) emission toward Sz 91 were conducted in 2010 
with the Submillimeter Array (SMA; Ho et al. 2004), which 
comprises eight 6 m antennas located atop Mauna Kea in 
Hawaii. The observations were performed in the compact 
configuration with baselines ranging from 6 to 70 m. Eor the 
continuum emission, higher resolution data were also obtained 
in 2010 in the very extended configuration with baselines up 
to ~500 m. The field of view of SMA was approximately 30'' 
in EWHM. The emission was detected using double side-band 
superconductor-insulator-superconductor receivers with a local 
oscillator frequency of Vlo = 340.755 GHz. Both the upper and 
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Figure 1. CLEANed maps of the 345 GHz continuum emission created by all UV data (a) and 50-600 kX (b). The contours begin at ±2cr with intervals of 2cr, 
where Icr =2.1 and 2.8 mJy beam“^ in (a) and (b), respectively. The red cross indicates the stellar position in each panel. The synthesized beam sizes in half power 
beamwidth (HPBW) are shown at the bottom-left corner in each panel: 174 x 170 with a position angle (P.A.) = — 9?8 and O'.'? x 0'.'3 with a PA. = 13?0 for (a) and 
(b), respectively. The blue box in the left panel indicates the area of the right panel. 

(A color version of this figure is available in the online journal.) 


lower sidebands data were used for the continuum, resulting 
in an 8 GHz bandwidth in total. The channel spacing of the 
CO(3-2) line was set to 0.8125 MHz. 

The amplitude and phase of the array system were cali- 
brated by observations of quasars, J1 604— 446, J1 626— 298, and 
J1517— 243, in a cycle with 12.5 min on a target and 6 or 8 min 
on a quasar. The response across the observed passbands was 
determined by 60 min observations of quasar, 3C 273. Absolute 
flux calibration was achieved by observing Titan at the begin- 
ning and end of each night. 

We created the continuum images by combining the very 
extended and compact configuration data, whereas the CO maps 
were created from the compact configuration data only. The 
UV data were edited and calibrated using MIR, an IDL-based 
software package. We used Astronomical Image Processing 
System for imaging procedures, including deconvolution by the 
CLEAN algorithm and restoration with a synthesized beam. 
Both continuum and line emission maps were created with 
natural weighting in the visibility plane, producing synthesized 
beam FWHMs of 174 x 170 at a position angle (P.A.) of —978 
and 476 x 178 at a P.A. of -470, respectively. 

2.2. High- re solution Polarized Intensity Imaging at 
the Ks Band with the Subaru Telescope 

^^-band (2.15 /xm) linear polarized intensity (PI) images of 
Sz 91 were obtained with the high-contrast imaging instrument 
(HiCIAO; Tamura et al. 2006) combined with dual-beam po- 
larimetry from the 8.2 m Subaru Telescope in 2012 May. The 
observations were conducted under the program SEEDS (Strate- 
gic Explorations of Exoplanets and Disks with Subaru; Tamura 
2009). The adaptive optics system (AO 188: Hayano et al. 2004) 
provided a limited diffraction and mostly stable stellar point 
spread function (PSF) with a FWHM of 0'.'06 in the Ks band. 
Polarization differential imaging (PDI) is a powerful technique 
used to reveal the structure of a dusty disk in very close proxim- 
ity to a star (e.g., Thalmann et al. 2010; Hashimoto et al. 2011, 
2012; Tanii et al. 2012; Muto et al. 2012; Kusakabe et al. 2012; 
Mayama et al. 2012). We employed the PDI mode, combined 
with the angular differential imaging mode, in which the field 


of view and the pixel scale were 10'' x 20" and 9.5 mas pixeU^ , 
respectively. Half-wave plates were placed at the four angular 
positions of 0°, 45°, 2275, and 6775 in sequence with a 30 s 
exposure per wave plate position. Image Reduction and Anal- 
ysis Facility (IRAF^^) software was used for all data reduction 
as following the methods of Hashimoto et al. (2011) and Tanii 
et al. (2012), and the Stokes Q and U parameter images were 
created. We calculated the PI as Q} and the polarization 
vector angle as 0.5 x arctan(I7 / Q) with the 3.4.0 version of the 
Common Astronomy Software Applications package. The final 
PI image was created from the smoothed Q and IJ images by 
a Gaussian function; thus, the effective resolution of the image 
became ~0'.'25. 

3. RESULTS 

3.1. 345 GHz Continuum Emission Map 

We created a 345 GHz continuum emission map of Sz 91 
with all UV data as shown in Figure 1(a). The map shows 
a strong continuum emission peak at the stellar position. The 
total flux density is measured to be 32.1 zb 3.6 mJy, which 
agrees well with the flux density of 34.5 zb 2.9 mJy, which 
is the previous 345 GHz measurement with the 12 m single 
dish APEX telescope obtained by Romero et al. (2012). The 
345 GHz emission is resolved by the 174 x I'.'O beam at a 
P.A. of —978. The beam-deconvolved size of the emission is 
measured to be ( 177 zb O'.' 1 ) x (O'.' 7 zb 0'.'2) with a P.A. of 1 697 0 zb 
370 from a two-dimensional Gaussian fitting, corresponding to 
(340 zb 20) X (140 zb 40) AU at a distance of 200 pc. This 
size is comparable to that of a typical protoplanetary disk. 
The 345 GHz emission is probably thermal dust emission from 
the disk around Sz 91 because the flux densities from 350 to 
1100/xm monotonically decrease with a spectral index of ~2 
(Figure 5), which is comparable to other T Tauri stars (Andrews 
& Williams 2005, 2007). 


IRAF is distributed by National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science Foundation. 
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Figure 2. CO(3-2) velocity channel maps with 0.35 km s ^ resolution obtained by SMA. The central LSR velocity of each panel is shown at the top in km s ^ . The 
contours start at ±2cr with intervals of 2cr, where Icr = 289 mJy beam“^. The ellipse at the bottom left in the top left panel shows the synthesized beam size of 
476 X 178 with a RA. of -4?0 in HPBW. 

(A color version of this figure is available in the online journal.) 


Moreover, we created a higher-resolution 345 GHz image 
using the UV data over the UV range from 50 to 600 kA. as 
shown in Figure 1(b). The emission peak at the 4.5a level can 
be found to the north of the star. The peak position is measured to 
be (ofj 2 ooo = 16^07”^! T.6, feooo = -39°03'47'.'8) and is shifted 
from the stellar position by 0'.'43 toward the direction at a P.A. 
of 1?4. The lack of a peak at the stellar position implies that 
the inner part of the dust disk is depleted or cleared as expected 
from the SED. Therefore, the continuum emission in Figure 1(b) 
probably originated from the innermost part of the dust disk. 
Assuming that the peak position in Figure 1(b) represents the 
inner edge of the dust disk, the inner radius of the Sz 9 1 transition 
disk is estimated to be 86 zb 25 AU. The total flux density of this 
compact component is measured to be 13.2 zb 3.1 mJy, which 
is 40% of the total flux density in Figure 1(a). 

3.2. CO( 3-2) Molecular Line Emission Maps 

Figure 2 shows velocity channel maps of the CO(3-2) line. 
The significant CO(3-2) emission is detected in the velocity 
range from 1.9 to 7.2 km s“^ in Vlsr- H is clear that the 
emission at Vlsr = 1.9-5. 8 km s“^ is concentrated at the 
stellar position. In addition, a velocity gradient appears to be 
present along the north- south direction in this velocity range: 
the emission is located mainly on the north side of the star at 
^LSR = 2.6-2.9 km s“^ and at the south from 4.0 to 5.1 km s“^. 
Notably, the central velocity range from 3.3 to 4.0 km s“^ 


includes the radial velocity of Sz 91 (Vlsr = 3.87 km s“\ 
converted from Vheiio = — 1 .57 km s“^ ; Melo 2003). In the range 
of 4.4-5.4 km s“^ , we detected spatially extended emission near 
the star (4.4 and 4.7 km s“^) and at the south side of the field of 
view (4. 7-5.4 km s“^). The origin of the extended components 
is most probably an ambient cloud because the systemic velocity 
and the velocity width of the main cloud condensation of Lupus 
III have been measured to be 4.1 and 1.2 km s“\ respectively 
(Kara et al. 1999). We also detected the extended emission in the 
6. 1-7.2 km s“^ channels. These components are shifted by ^3" 
toward the south-west direction from the star and also possibly 
originate from the spatially extended ambient cloud. Hereafter, 
we focus on compact components toward the star that should 
have originated from the gas disk; we do not discuss the extended 
components in this paper. 

To clearly see the gas disk emission and its velocity gradi- 
ent, we created total integrated intensity (1. 5-6.1 km s“^) and 
first moment maps of the CO(3-2) line, as shown in Figure 3. 
From these maps, we successfully detected the centrally con- 
centrated emission and the velocity gradient roughly along the 
north-south direction, suggestive of the rotating gas disk. The 
total integrated intensity of CO(3-2) in the area above the 3a 
noise level is 8.55 Jy km s“^ . The emission peak is shifted from 
the stellar position by toward the north direction, possibly 
because the red-shifted emission is partially resolved-out owing 
to the contamination by the ambient cloud. 
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Figure 3. Total integrated intensity map of the CO(3-2) emission over a Vlsr range of 1 .5-6. 1 km s“^ (contours) superimposed on the intensity-weighted first moment 
map (color). The contours start at ±3cr with intervals of 3cr, where Icr = 450 mJy km s“^ beam“^ . The color bar on the right-hand side of the panel shows the centroid 
velocity in m s“^ . The red cross indicates the stellar position. The synthesized beam size in HPBW is shown at the bottom left corner: 476 x 178 with RA. = — 4?0. 
(A color version of this figure is available in the online journal.) 
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Figure 4. (a) Polarized intensity (PI) image obtained with the Subaru Telescope in the Kg band (color) where the 345 GHz continuum emission map (contour) in 
Figure 1(b) is superposed. The stellar position is shown by the red cross. The best-fit ellipse is indicated by the yellow dashed line. The inner area of (y.'25 diameter 
(filled circle) is photometrically unreliable owing to the point spread function (PSF) subtraction process and is masked, (b) Polarization vector angle map superimposed 
on the PI image in the left panel. The vector directions indicate angles of polarization. The vector’s lengths are arbitrary. 

(A color version of this figure is available in the online journal.) 


3.3. Polarized Intensity Image at the Kg Band 

The high-resolution PI and polarization vector angle maps 
at the .^^-band are shown in Figures 4(a) and (b), respectively. 
The presented images are smoothed and the effective resolution 
of the images is 25. A crescent-like emission region was 
detected around Sz 91, which is elongated from the south to 
the north via the west side of the star. The polarization angles 
are nearly perpendicular to the radial directions from the central 
star, indicating that the ^^-band emission probably originated 
from the scattered light at the inner part of the dust disk. The 
crescent-like emission suggests that this part is the near side of 
the disk if we assume that forward scattering is dominant, as 
is the case in Mie scattering. Substantial emission dips appear 
at the north and south sides of the software mask, possibly due 
to the inner hole structure of the dust disk. The PSF subtraction 


process is a primary factor to cause an artificial systematic error 
near the star, and the error becomes a systematic emission via 
the square root operation of stokes Q and U. The symmetric 
distribution of the dips indicates the existence of an inner hole 
in the dust disk. 

We also detected a bright region near the western edge of 
the software mask, i.e., a photometrically reliable emission. 
However, the origin of this emission remains unclear. Because 
the PSF shows a substantially symmetric distribution, it is 
unlikely that the bright emission would appear at the only one- 
side. Because it extends to the edge of the software mask, higher 
resolution and higher contrast observations close to the star are 
required to determine the origin of the bright emission. 

The existence of the inner hole structure expected from the 
PI image strongly supports the results of the submillimeter 
images. The higher-resolution 345 GHz image is overplotted in 
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Table 1 

Parameters of the Best-fit Ellipse in Figure 4 


Parameter 

Best-fit Value 

Error 

AR.A.^ (") 

0.03 

0.04 

ADecL^ CO 

0.03 

0.02 

Major axis ('0 

0.33 

0.02 

Minor axis (") 

0.25 

0.08 

Position angle (°) 

17.5 

17.7 


Note. ^ Offset from the stellar position. 


Figure 4(a). It is clear that the PI is distributed at the innermost 
part of the submillimeter emission. For example, along P.A. = 
0°, the 345 GHz emission peaks at (X.'45 from the star whereas 
the PI peaks at O'.' 32. These facts suggest the existence of the hole 
structure at the inner part of the disk and that the PI originates 
from the innermost part of the transition disk, i.e., the inner 
edge. 

To quantify the shape of the inner edge of the dust disk, 
we fit an ellipse to the crescent-like emission. Table 1 lists the 
parameters of the best-fit ellipse, and the fitted ellipse is shown 
in Figure 4(a) by the yellow dotted line. The center position 
of the ellipse coincides with the stellar position. The P.A. of 
the ellipse, i.e., the direction of the disk major axis, is along at 
17? 5 zb 17? 7, which is consistent with the P.A. of the emission 
dips in the crescent-like emission. The PI at the disk minor axis 
(P.A. 120°) appears to be lower than that at the major axis. 

This result is probably due to lower polarization degree along 
the minor axis where the scattering angle at the disk surface is 
deviated from 90° (McCabe et al. 2002). 

4. DISCUSSION 

4.1. Spatially Resolved Disk Structure 

Both the high-resolution submillimeter and NIR images 
clearly resolved the disk structure around Sz 91. In this section, 
we estimate the inner radius, outer radius, and inclination angle 
of the disk directly from our images. These parameters are 
adopted in the model calculations in Sections 4.2 and 4.4. 

The inner radius of the dust disk can be estimated from the 
large inner hole revealed in the ^^-band PI image. Assuming that 
the -band emission is the scattered light at the inner edge of 
the disk, the best-fit ellipse in the /T^-band image shows the inner 
radius to be 65 zb 4 AU. There are two possibilities to explain 
the scattered light. One is forward scattering at the surface of the 
inner part of the disk, and the other is scattering at the inner wall 
of the disk such as that in the case of LkCa 15 (Thalmann et al. 
2010). We can not judge which case is more likely to explain the 
NIR emission from the Sz 91 disk. However, in either case, it 
probably originates from the innermost part of the disk, and our 
estimation of the inner radius seems to be reasonable. The high- 
resolution 345 GHz image also provides a rough estimate of the 
disk inner radius. Although the 345 GHz image in Figure 1(b) 
does not clearly show the presence of the inner hole, it is clear 
that the 345 GHz flux density becomes weak toward the central 
star, which is suggestive of dust depletion near the star. The 
angular distance between the 345 GHz emission peak and the 
central star can be interpreted as the disk inner radius of (y.'43 
(=86 AU), which is roughly consistent with the above estimate 
of 65 AU. In the following discussion, we adopt the disk radius 
Rin of65 AU. 

The outer radius of the dust disk can be estimated from the 
beam-deconvolved size of the full-UV 345 GHz map shown 


in Figure 1(a). Assuming a geometrically thin disk, the beam- 
deconvolved size of (1'.'7 zb (y.'l) X ((y.'7 zb (y.'2) corresponds to 
the outer radius Rout of 170 zb 20 AU at the distance of 200 pc. 
The gas disk is also resolved in CO(3-2) (Figure 3). The beam- 
deconvolved size is measured to be (3'.'1 zb 0'.'3) x (1'.'6 zb 0'.'7) 
at a P.A. of 15 zb 11° from a two-dimensional Gaussian fitting, 
corresponding to Rout = 310 zb 60 AU. However, the red-shifted 
gas is probably contaminated by the ambient cloud gas, and the 
center position of the beam-deconvolved disk is shifted by 0'.'6 
from the star to the north. Thus, we consider the largest distance 
of 2'.' 1 from the stellar position as the outer radius of the gas 
disk: Rout = 420 AU. Because the actual intensity distribution 
of each emission is expected to have a power-law form, these 
outer radii derived from the Gaussian fitting should be regarded 
as lower limits. 

The outer radius of the gas disk is much larger than that 
of the dust disk. The discrepancy in radial extension between 
the dust and gas emissions is frequently observed for other 
protoplanetary disks, which can be explained by an exponential 
decrement in the surface density (e.g., Hughes et al. 2008; Panic 
et al. 2009; Isella et al. 2010) or by the radial drift of large 
grains (Andrews et al. 2012). In addition, the possible north- 
south asymmetry in the dust continuum emission, suggested 
in Figure 1(b), could have resulted from the azimuthal drift of 
large grains (Bimstiel et al. 2013; van der Marel et al. 2013) 
associated perhaps with a perturbing body orbiting within the 
central hole of this disk. 

The disk inclination angle can be derived from the ratio of the 
major to minor axis lengths of the best-fit ellipse in Figure 4(a): 
40 zb 15° from face-on. Moreover, the beam-deconvolved size 
of the 345 GHz continuum emission in Figure 1 (a) provides an 
inclination angle of 66 zb 5°. The latter value of 66°, however, is 
possibly influenced by the disk flaring because the distribution 
of the 345 GHz emission is biased to the outer part of the disk. 
Thus, we adopt the R^-band value of 40° as the disk inclination 
angle in this study. In fact, for the CO(3-2) spectrum, the model 
calculation in Section 4 prefers the smaller value of 40°. 

4.2. Disk Parameters Derived from SED Fitting 

To deduce other parameters of the dust disk, we performed 
least- square fitting to the SED on the basis of the power-law 
disk model. Figure 5 shows the SED of Sz 91 that includes our 
345 GHz flux density and a recent data set by WISE. As Romero 
et al. (2012) reported, the SED is characterized by a large dip 
around 20 /xm, a sharp rise from 20 to 70 /xm, and a significant 
(sub-)millimeter emission. 

We introduced the two components of a cold disk and a small 
amount of hot dust inside the disk as follows. First, we applied 
the model for the cold disk, which comprised the usual power- 
law disk and blackbody stellar emissions, to all of the SED data 
except that at 10-30 /xm, because the model could not reproduce 
all the data points including the three points at 10-30 /xm. The 
power-law disk model has surface density and temperature radial 
distributions of S(r) and T(r), respectively, in a power-law form 
(Kitamura et al. 2002; Tsukagoshi et al. 2011): 

= ( 1 ) 

and 

(2) 
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Table 3 

Best- fit Parameters in the SED Fitting 
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X (iLim) 

Figure 5. Spectral energy distribution of Sz 91. The filled circles indicate 
observed flux densities. All data except the 345 GHz flux density were compiled 
from previous studies, including the NOMAD catalog (Zacharias et al. 2005), 
the 2MASS point source catalog (Cutri et al. 2003), the Spitzer IRAC and 
MIPS photometry (Evans et al. 2009), the WISE all-sky data release (Cutri et al. 
2012), the AKARI FIS all-sky survey point source catalog (Yamamura et al. 
2010), the IRAS point source reject catalog (Infrared Astronomical Satellite 
(IRAS) Catalogs 2007), the AzTEC/ASTE 1100/xm photometry (R. Kawabe 
et al. 2014, in preparation), and the Herschel PACS/SPIRE FIR flux densities 
which we derived from the HSA Science archive data. The best-fit SED and 
the contribution of the hot component in the optically thin case are shown by 
the black and red lines, respectively, and by the gray and orange lines in the 
optically thick case. The blue dashed line indicates the stellar contribution. 

(A color version of this figure is available in the online journal.) 


Parameter 

Value 

Stellar parameters 

{Rq) 

1.29 ±0.03 

n(K) 

4148 ± 55 

(^o) 

0.49 ± 0.02 

Cold outer disk 

fin (K) 

32.5 ± 3.9 

Zin^ (g cm“2) 

0.67 ± 0.03 

h 

0.5 ±0.1 

Mdisk" (10-3 Mq) 

2.4 ± 0.8 


Note. ^ The surface density and the disk mass 
are shown in the gas+dust density and mass 
by assuming the gas-to-dust mass ratio of 100. 


component is supported by the sign of mass accretion derived 
from the Ha emission line (Romero et al. 2012). Because no 
information was available on the structure of this component, 
we simply assumed a gray body with a temperature of Tc, a 
column density of Ec, and a solid angle of Qc- Us flux density 
Sx at wavelength X is written by 


Table 2 

Fixed Parameters in the SED Fitting 


Parameter 

Fixed Value 

M, (Mo) 

0.49^ 

Av (mag) 

2.0^ 

Rm (AU) 

65 

^out (AU) 

170 

i (deg) 

40 

P 

1.5 

q 

0.5 


Note. ^ Hughes et al. (1994) 


where r is the radial distance, Ein and Tin are the surface density 
and temperature at the inner radius, Rin, and p and q are power- 
law indexes. The gas-to-dust ratio was assumed to be 100 and 
the extinction of the stellar light because of the interstellar dust 
was corrected by the Ay value toward the star and the dust mass 
absorption coefficient (Figure 1 of Adams et al. 1988). Table 2 
lists the fixed parameters in the SED fitting. The stellar mass, 
the effective temperature of the star, and the visual extinction 
toward the star are from Hughes et al. (1994). The inner and 
outer radii and inclination angle of the disk were determined in 
Section 4.1. We here adopted the values of 1.5 and 0.5 as p and 
q, respectively, which are the same as those of the minimum 
mass solar nebular (Hayashi 1981). These values were selected 
because the SED fitting is known to be insensitive to the power- 
law index of p, and it is difficult to resolve the parameters for the 
temperature profile {T and q) from only the longer wavelength 
data (>30 /xm). Notably, the lower limit of the disk temperature 
was set to be 10 K, which is the typical temperature of the Lupus 
III cloud (Vilas-Boas et al. 2000). The stellar radius, Tin, Ein, 
and the power-law index of the dust mass opacity coefficient ^ 
{Ky = 0.1 X (v/10^^ Hz)^ cm^ g~^: Beckwith et al. 1990) were 
treated as the free parameters in the SED fitting. 

After the power-law disk model fitting, we introduced the 
additional hot component inside the disk to reproduce the 
observed flux densities at ^^20 /xm. The presence of the inner hot 




2hP 1 


(1 - X 


( 3 ) 


where h is the Plank constant, c is the speed of light, and Xb is 
the Boltzmann constant. We applied the same dust mass opacity 
coefficient, /Cy, as that in the cold disk. The parameters of Tc, Ec, 
and Qc were treated as free. In the fitting, the upper limit of Qc 
was set to be ~6.0 x 10“^^ str, corresponding to the solid angle 
of the inner hole seen in the NIR image. 

The best-fit model SED of the cold disk reproduced the longer 
wavelength data (> 30 /xm) effectively, as shown in Figure 5. The 
best-fit parameters of the cold disk are summarized in Table 3. 
The P value of 0.5 zb 0.1 is significantly smaller than that in 
the diffuse interstellar medium of ^2 (Draine & Lee 1984), 
suggesting that the dust growth occurs in the disk (Miyake 
& Nakagawa 1993). By adopting this ^ and Ky introduced 
by Beckwith et al. (1990), who assumed a 100:1 mass ratio 
between gas and dust, the disk mass is derived to be (2.4 zb 
0.8) X 10“^ Mq. This value is significantly higher than those of 
other class III sources in nearby star forming regions (Andrews 
& Williams 2005, 2007); however, the value is lower than those 
of most of the transition disks studied thus far (Andrews et al. 
2011), even if we consider the difference in Ky by a factor of 
~L7 at V = 340 GHz between this study and the previous 
studies in which ^ is set to be 1 . 

To verify the fitting result, we also created a continuum image 
for the best-fit model and compared it with the observations, as 
shown in Figure 6. The model image essentially agrees with the 
observations, but there remains a difference of at most 4a; the 
negative residual extends in the east-west direction, which may 
be due to the asymmetry of the disk. 

We next attempted to reproduce the observed flux densities at 
10-30 /xm by adding the contribution from the hot component. 
The best-fit parameters of Tc, Ec, and Qc to reproduce all of the 
SED data were searched in the reduced maps by manually 
changing the initial values of Tc, Ec, and Qc with 6000 runs. 

Erom the calculations, we determined that there are two 
distinct regions in the Ec-Qc plane where the reduced x^ 
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R.A. (J2000) 

Figure 6. Comparison of the observed data with the best-fit model. The 345 GHz continuum image, model image, and residual image (data-model) are shown from 
left to right. The contour starts at ±2cr, where Icr =2.1 mJy beam“^ The solid and dashed lines indicate positive and negative, respectively. 


takes its local minimum values, which indicates a strong 
coupling between Ec and Qc- One is the region in which 
Ec < 5 X 10“^ g cm“^ and Qc > 4 x 10“^^ str, and the 
emission of the hot component is substantially optically thin. 
In this region, the best-fit typically converges to be 186 K 
and the best-fit Ec is inversely proportional to Qc. The product 
of Ec X Qc is a constant value, which provides the best-fit mass 
of 3 X 10“^ Mq for the hot component. The other is the region 
in which Ec > 30 g cm“^ and the emission is substantially 
optically thick. In this region, the best-fit is typically 172 K 
and Ec is independent of Qc. For the hot component, we can 
obtain the best-fit Qc = 1 x 10“^^ str and the lower limit mass 
of 6 X 10“^ Mq. Although the reduced values are slightly 
lower in the former case (Ax ^ = 0.3), the difference between the 
two regions is not significant. We therefore conclude that there 
are two possible origins of the hot component: the optically 
thin gray body emission and the optically thick black-body 
emission. However, the optically thin condition is the unlikely 
case because the mass of the hot component is too small given 
the mass accretion rate of 10“^^ Mq yr“^ ; all the hot component 
would disappear in only 30 yr. 

Notably, the contribution of the hot component is mainly 
restricted by only three data points at MIR, from 12 to 25 /xm; 
the upper limit of the temperature is determined by the slight 
excess emission at 12 /xm and the lower limit is limited by 
the decrement between 22 and 24 /xm. A refined model will be 
required when more data at MIR and FIR are obtained by further 
observations. 

4.3. Origin of the Hot Component in the Inner Hole of the Disk 

The SED model analysis indicates that the presence of a hot 
component inside the dust disk with a temperature of ~ 180 K, 
which is responsible for the SED peak at "^20 /xm. There are 
two possible origins for such a hot component: a localized self- 
luminous emitting body (i.e., circumplanetary disk) or an inner 
warm structure of the disk (Wolf & D’Angelo 2005). In this 
section, we estimate the size and mass of the hot component for 
both cases. 

In the case of the circumplanetary disk, we can estimate the 
radius from Qc because it is usually optically thick (D’Angelo 
et al. 2003). Assuming the circumplanetary disk is parallel to 
the parent circumstellar disk, the solid angle of 1 x 10 str 
corresponds to the radius of 0.3 AU, or 64 Rq . On the other hand. 


the Hill radius of a putative planet around Sz 91 is expressed by 

j?Hiii = 1.1 X ^ [AU], (4) 

lOAuVlM,/ 

where m and a are the mass and orbital radius of the planet, 
respectively. The radius of the circumplanetary disk inferred 
from the SED fitting is significantly smaller than Rhiii at 
3-65 AU from the star, which is consistent with the theoretical 
expectation for a circumplanetary disk (e.g., Tanigawa et al. 
2012). However, our angular differential imaging observation 
with the Subaru Telescope could not examine the presence of 
such a companion planet due to a low rotation angle of ~16°. 

For the second possibility, because the hot component could 
be fitted with a single temperature, it must be confined to 
a narrow width in radius. If we extrapolate the temperature 
distribution of the outer disk determined in the SED fitting, 
the temperature of the hot component in the optically thick 
case (172 K) corresponds to a radius of 2.3 AU. The best-fit 
Q^ = 1 X 10“^^ str corresponds to a ring width of 0.01 AU, 
which is significantly narrow with respect to the ring radius. 
In contrast, the optically thin condition is unlikely because the 
solid angle is >4 x 10“^^ str which corresponds to a ring width 
of >53 AU. Therefore, we conclude that the other implication 
of the hot component is the optically thick ring at 2.3 AU, 
whose total mass is at least 6 x 10“^ Mq. Such an example of 
the optically thick ring around the transition disk has also been 
reported in RX J1633. 9-2442 (Cieza et al. 2012a). The inner 
structure of Sz 91 may be similar to that of RX J 1633.9-2442, 
whereas Sz 9 1 exhibits a larger inner hole and a lower disk mass. 

4.4. Velocity Structure of the Gas Disk: Model Calculation and 
Comparison with the CO( 3-2 ) Profile 

The CO(3-2) image shown in Figure 3 suggests the presence 
of a rotating gas disk around the star. To reveal the disk rotation 
in detail, we calculated model spectra of CO(3-2) with a simple 
power-law disk model according to that reported by Kitamura 
et al. (1993), and we compared the results with the observed 
CO(3-2) profile. 

The following parameters of the model disk were estimated 
from the observed images and the SED fitting: the inner and 
outer disk radii of 65 and 420 AU, respectively; the temperature 
distribution of T{r) = 32.5 x {r / K; and the surface 
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Figure 7. CO(3-2) spectrum integrated over a 5" x 5" area centered at the stellar 
position (black line). The red line shows the calculated profile by a disk model 
with an inner radius of 65 AU and an inclination angle of 40° from face-on (see 
the text for the detailed description of the model). The blue line indicates the 
case of an inner radius of 0 AU. The green line indicates the model profile of 
the disk with a 65 AU inner radius and an inclination of 66° . The CO emission 
in the velocity range from 4 to 7 km s“^ in Vlsr is probably contaminated by 
the ambient cloud emission, as shown by the gray line at the top of the panel. 
(A color version of this figure is available in the online journal.) 


density distribution of E(r) = 0.67 x g cm~^. The 

hydrostatic equilibrium is assumed along the vertical direction 
and the density distribution, p(r, z), is therefore expressed by 


p(r, z) = p(r, 0) exp 




2n 


(5) 


where H(r) is the scale height given by 


H(r) = 


r^k^T{r) 
GM^jjimu ' 


( 6 ) 


Here, G is the gravitational constant, is the stellar mass, [i 
is mean molecular weight, and mn is the mass of H atom. The 
density at the midplane is expressed by 


P(r, 0) = 


^r) 

V^H(r) ' 


(7) 


The gas motion is assumed to be Kepler rotation whose velocity 
field is written by 



The fractional abundance of CO with respect to H 2 is assumed 
to be 9 X 10“^, which corresponds to a typical interstellar value 
(e.g., Schloerb & Snell 1984; Irvine et al. 1985). Because the 
density of the disk is at least 10^ cm“^ at the midplane which 
is significantly higher than the critical density of CO(3-2) 
(~10^ cm“^), we assume local thermodynamic equilibrium. 
Notably, the hot component in the dust hole was not included 
because its contribution was negligible (^1 x 10“^ Jy). 

Figure 7 shows the calculated CO spectra superimposed 
on the observed CO spectrum integrated over a 5'' x 5'^ box 
centered at the stellar position. We noted that the peak intensity 
of the calculated CO spectrum differed from the observed 


value. However, the discrepancy is not significant because 
the disk emission around the stellar LSR velocity is probably 
contaminated by the ambient cloud emission. Therefore, we 
focused on the emission at the blue-shifted side (Vlsr < 
2.9 km s“^). The model profile in the i = 40° case agrees with 
the observed data, confirming the validity of the disk parameters 
derived from the dust disk. Notably, the inclination angle of 66° 
determined by the 345 GHz image did not fit well the CO line 
shape, indicating that the inclination of the dust disk estimated 
from the NIR scattered light is more plausible. 

The presence of the gas inside the inner edge of the dust 
disk could not be confirmed from our data set. Although we 
calculated the model spectrum for Rin = 0 AU as shown 
in Figure 7, the difference is within the Icr uncertainty. The 
presence of the gas disk in the inner hole is supported by the 
fact that Sz 91 shows the mass accretion onto the star, and thus, 
higher spatial resolution imaging is required to reveal the inner 
structure of the gas disk. 


5. SUMMARY 

We present the results of the aperture synthesis 345 GHz 
continuum and CO(3-2) line emission observations with SMA 
toward a transition disk object in Lupus, Sz 91. Furthermore, a 
high-resolution image of the PI at the Ks band obtained with the 
Subaru Telescope is also presented. The transition disk around 
Sz 91 has been directly resolved and imaged in this study. The 
disk parameters are derived and the structure and evolutional 
phase of the Sz 91 transition disk are discussed. The main results 
of our observations are summarized in the following points: 

1. Our high-resolution imaging revealed a dust disk around 
Sz 91 with inner and outer radii of 65 and 170 AU, 
respectively, and an inclination angle of 40° . Furthermore, 
the Kepler rotating gas disk with a radius of 420 AU was 
imaged in the CO(3-2) line. 

2. Model analysis of the SED of Sz 91 was performed using 
a simple power-law disk model. We determined that the 
observed SED can be reproduced well by the combination 
of a cold disk and a hot component in the inner hole of 
the disk. The total H 2 mass of the cold disk is estimated to 
be 2.4 X 10 ^ Mq if the canonical gas-to-dust mass ratio 
of 100 is adopted. The disk mass is significantly higher 
than those of other class III sources in nearby star forming 
regions; however, the disk mass is one of the lowest masses 
among the currently known transition disks. 

3. We determined that the hot component can be expressed by 
a single temperature gray body of ~ 1 80 K. Although the hot 
component could not be resolved by our observations, its 
origin is either a localized self-luminous emitting body (i.e., 
a Jovian mass protoplanet with a circumplanetary disk) or 
an optically thick ring in the inner hole of the disk at 2.3 AU. 

4. Our results confirm the previous results such that the disk 
structure of Sz 91 is consistent with that of an ongoing giant 
planet forming disk. In particular, the relatively large inner 
hole and lower disk mass indicate that the transition disk 
of Sz 91 is probably in a stage of nearly completing planet 
formation. Sz 91 will be a crucial target for investigating 
the evolution of transition disks and the planetary formation 
process. In the near future, our proposed study with the 
Atacama Large Millimeter/ submillimeter Array (ALMA) 
will provide a new insight into the planet formation process. 
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